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Abstract. We consider a discrete-time quantum walk Wt, K at time ( on a graph with joined half lines J K , which 
is composed of k half lines with the same origin. Our analysis is based on a reduction of the walk on a half line. 
The idea plays an important role to analyze the walks on some class of graphs with symmetric initial states. In this 
paper, we introduce a quantum walk with an enlarged basis and show that Wt, K can be reduced to the walk on a 
half line even if the initial state is asymmetric. For Wt, K , we obtain two types of limit theorems. The first one is 
an asymptotic behavior of Wt,n which corresponds to localization. For some conditions, we find that the asymptotic 
behavior oscillates. The second one is the weak convergence theorem for Wt, K - On each half line, Wt,n converges to 
a density function like the case of the one-dimensional lattice with a scaling order of t. The results contain the cases 
of quantum walks starting from the general initial state on a half line with the general coin and homogeneous trees 
with the Grover coin. Q 

1 Introduction 

Random walks have a very important role in various fields, such as physical systems, mathematical modeling 
and computer algorithms. In 1990s, quantum walks arise as a quantum counterpart of random walks [TH3]. 
They are defined by unitary evolutions of probability amplitudes, whereas random walks are obtained by 
evolutions of probabilities by transition matrices. Discrete-time quantum walks are introduced by Refs. [H2]- 
In recent years, quantum walks have been well developed in fields of quantum algorithms, for example [3HS]- 
On the other hand, studies of the walks from the mathematical point of view also arise. Especially, as 
a limiting behaver, localization appears in quantum cases [7r412j. Furthermore the quantum walk has a 
quadratically faster scaling order than the random walk in the weak convergence [T3r417| . Cantero et al. 
introduced an analysis using the CMV matrix [TTJ[T2]. This method is very useful to consider localization. 
To analyze the quantum walk, we use the generating function. By using the generating function, we can 
compute not only localization but also the weak convergence of the walk. A reduction technique [2DH^] . 
which reduces the walk to a one-dimensional quantum walk, is very important to apply a path counting 
method [13, 14 , 23J which gives an explicit expression for the generating function. To treat the quantum walk 
with asymmetric initial states, we introduce a quantum walk with enlarged bases. 

Our main results are two limit theorems for the quantum walk Wt. K on a graph with joined half lines 
with arbitrary initial state starting from the origin. In case of k = 1, Wt,i corresponds to a quantum walk 
on a half line with the general coin. Furthermore, by considering the reduction of the walks, the two limit 
theorems can be adopted to quantum walks on homogeneous trees and semi-homogeneous trees with the 
Grover coin operator. One of two our main results is the explicit expression for the limit probability of 
Wt, K - It is corresponding to localization which is defined that there exists a vertex of the graph x such that 
limsupt^^ P(Wt, K — x) > 0. We find that, for some conditions, the asymptotic behavior oscillates. Same 
as other results on quantum walks [7H2], localization has an exponential decay for position x on each half 
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line. Another main result is the weak convergence of Wt }K - On each half line, Wt, K has a scaling order t. 
Moreover the limit measure has a typical density function which appears on other quantum walks [8, 13HT7] . 

For related works, Chisaki et al. jSj obtained the same type of limit theorems for a quantum walk on 
homogeneous trees with two special initial states. This result induces limit theorems for a quantum walk on 
a half line with a special coin operator. Konno and Segawa |18] showed localization of quantum walks on a 
half line by using the spectral analysis of the corresponding CMV matrices. 

The remainder of the present paper is organized as follows. In Section 2, we give definitions of discrete- 
time quantum walks treated in this paper. Section 3 presents our results. Section 4 gives proofs of our 
main theorems. In Subsection 4.1, we introduce a quantum walk with an enlarged basis and reduce Wt, K 
to the walk on a half line. Subsection 4.2 presents a proof of Theorem 1 based on the generating function. 
Subsection 4.3 is devoted to a proof of Theorem 2 using the Fourier transform of the generating function. 
In Appendix, we compute the generating function. 



2 Discrete-time quantum walks 

This section gives the definition of the quantum walk on undirected connected graph G. Let V(G) be a 
set of all vertices in G and E(G) be a set of all edges in G. Here we define E X (G) C E(G) as a set of 
all edges which connect the vertex x £ V(G). Now we take a Hilbert space spanned by an orthonormal 
basis {\x); x £ V(G)} as a position space Jt? p and a Hilbert space generated by an orthonormal basis 
{\l); I £ E X (G)} for x £ V(G) as a local coin space J^ Cx - A discrete-time quantum walk on G is defined 
on a Hilbert space ffl spanned by an orthonormal basis {\x,l); x £ V(G), / £ E X (G)}. Note that if we 
take G as a regular graph, J$? can be written as Jrf? = Jtfp ® JO^ for any x £ V(G). On the space Jtf, the 
evolution operator U is given by U = SF, where S : Jf? —> is a shift operator and F : iff — > iff is a coin 
operator. Here we define F — J2 x ev(G) \ x )( x \ ® ^ as a coin operator and C x : — » for x £ V(G) as 
a local coin operator. If the graph is regular and the local coin operator is all the same, we can rewrite the 
coin operator as F = I p <g> C, where I p is the identity operator on J$f p . As typical local coin operators, the 
Hadamard operator H and the Grover operator Gd are often used, where H and Gd {d > 2) are defined by 



Gd = 



ad b d 
bd a d 



bd b d 



b d 



In this paper we define a\ = 1, b\ = 2 and G\ 
x is described as 



1. From the construction, the state at time t and position 



*t(s) = J2 ^t{x,l)\x,l), 



(2.1) 



where «t(x, I) £ C is the amplitude of the base \x, I) at time t and C is the set of all complex numbers. The 
probability of the state is given by a square norm of ^t(x), i.e., ||^ r t (s)|| 2 = J2ieE x (G) I^^OI 2 - We only 
consider the initial state starting from the origin "o" with the state ^o{o) such that 1 1 M/q (c) 1 1 = 1- 



2.1 Quantum walk on a graph with joined half lines 

This subsection gives the definition of a graph with joined half lines J K and the quantum walk Wt. K on JJ K . 
Let K K = {0, 1, . . . , K - 1} and 7L r = {h r (l), h r (2), . . .} for r £ K K , we define V(J K ) = {0} U {U^k^-}. A 
vertex hi(x) connects hj(y) if and only if \x — y\ = 1 with i = j, and the origin connects h r (l) for any r 
(see Fig. [T] (a) for example). 

The quantum walk on J K is defined on Jf?^ K > which is a Hilbert space spanned by an orthonormal basis 
{|0, 1); I £ {eo, ei, . . . , e«-i}} U {\x, I); x £ V(I K ) \ {0}, I £ {Up, Down}}. Throughout this paper, we put 
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the base \e r ) as T [5 • • • 1 • • • 0], where T is the transposed operator. We define a local coin operator C as 



C = 



a b 
c d 



G U(2) with abed ^ 0, 



where [/ (d) is the set of d x d unitary matrices. The coin operator Fj is given by 

Fj = \0){0\®G K + J2 \x){x\®C, 

where G K is the Grover operator. The shift operator Sj is given by 



(2.2) 



(2.3) 



Sj\0,l) = 1^(1), Dowi) , / = e r , 

|0,e r ), / = C/p, 

|/i r (2), Down), I = Down, 

\h r (x-l),Up), l = Up, x>2 
\h r (x + 1), Down), I = Down, ~ 



Sj\h r (l),l) = 
Sj\h r (x),l) = 



Then the evolution operator of the walk Uj is obtained by Uj = SjFj. An expression of Wt,3 using weights 
is shown in Fig. [T] (a) , where 



Pi 1 



P«2 



P«3 



P = 



a b 
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a b 

a b 
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Note that W^i is a quantum walk on a half line with a reflecting wall, 
(a) (b) 





Figure 1: (a) Quantum walk on J3, (b) Quantum walk on T3 
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2.2 Quantum walk on homogeneous trees 

We define a homogeneous tree T K and a quantum walk Vt, K on T K 
be the set of generators subjected to the relation a 2 — e for j G K 



Fix k > 2, let S = {<7o, <j\,..., ct k _i} 
where the empty word e is the unit of 



this group. Then we put V(T K ) = {e} U {<Ji n . . . Oi 2 Oi 1 : n > 1, tr^. € E, ij+x ^ ij for j = 1, 2, . . . , n — 1}. 
Here vertices g and /i are connected if and only if ghr 1 G £. On this graph, JfCp is generated by an 
orthonormal basis {\g); g G V(T K )} and Mc" K is associated with an orthonormal basis { | <r j ) ; <jj G £}. We 
choose G K as the local coin operator, then the coin operator Ft and the shift operator St are defined as 
follows: for a G £ 



F T = \e){e\®cG K 
St\9,o) = \<rg,<r), 



E 

geV(T K )\{e} 



where we put |<j r ) as T [cT • -010- • -0] and c G C with |c| = 1. The phase c works as a defect on the origin, 
which is an extension of our model in [8]. An expression of Vj^ using weights is shown in Fig. Q] (b), where 



_I 2 2 

3 3 3 







, Q 3 = 





2 _ 1 

3 3 





i?3 = 








and P 3 = 5P 3 , Q 3 = cQ 3 , R 3 = SR 3 . 

In the case of the one point initial state on the origin, Vt >K can be reduced to the equivalent walk on 
JJ K even if the initial state is not symmetric. To explain it, we define subgraph G T K as V(T^) = 
{(Ji n ■ ■ ■ (Ti 2 ai 1 : n > 1, o\ . G S, cr^ — o~ r , ^ ij for j = 1, 2, . . . , n — 1}. They are subtrees whose roots 
are the children of the root of T K . Now we consider the following new basis, for x > 1, 



i^^p>^ = v 1 E E i*«*>. 

\g\ = x 



\g,°i) 



g 6 VCT^') <Tj:|<Tj3|=a:+l 

\g\ = x 



The new space Jif^"^ spanned by a basis {\e,l) : Z G £} U {\x, l) (Tr : r G K K ,a; G Z + ,Z G {[/p, Down}} is 
isomorphic to Jf?^^ under the following one-to-one correspondence 



|x, l)cr r «-» |/ir(ic), for / G {[/p, Down}, x > 1, 
|e,(T r ) -H- |0,e r ), 

where Z + = {1,2,...}. Then the direct computation gives the following lemma. 



(2.4) 



Lemma 1 (Homogeneous tree) The subspace of Jf^ K ^ is invariant under the action of the time evolu 
tion ofVt, K . In particular, when we take the bijection from Jf?^^ to Ji?^^ given by Eq. \2.4\ ), the walk 
equivalent to Wt, K with the following coin operator 



is 



ieV(l«)\(o} 



When we consider c~ 1 Fj , the above equation becomes a special case of Eq. (|2.3[) . since |c| = 1. 

Similar to Vt >K , we can define a quantum walk Vt, K ' tK on a semi-homogeneous tree T K i >K , which is a k- 
regular tree except the origin whose degree is k' > 2, with the local coin operator cG K ' at the origin and G K 



Vk - Tb K 



\Jk— lb K 
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otherwise. Then we can reduce it to Wt tK ' with the coin operator Fj K given by 



FY = |0> <0| OcG K /+ E \x){x\® 

i£V(l s /)\{0} 



y/n-lb K -a K 



The infinite binary tree is a special case for this graph (k = 3, k' = 2). 

3 Main results 

In our main theorems, we give explicit formulae with respect to each half line in J K . Let "ft (x) be the state 
of the quantum walk Wt, re at time t and position x. For x € Z + U {0}, we introduce random variables Xt >r as 
P(X t ,r = 0) = |a t (0,e r )i 2 and P(X t , r = 0) = \\V t (h r (x))\\ 2 . Remark that P(W t , K = h r (x)) = P{X t>r = x) 



for x > 1, P(W tlK = 0) = E,- £Kk = 0) and E, 



cGZ+U{0} 



P(X 



1. 



In order to describe the limit theorems for Wt, K , we first introduce several parameters. 

<p = arg (c), 

K± = \l±c\\ 

K x = (l- c )(l + c), 

where a is the complex conjugation of a E C. Next, we denote the following notations to state Theorem [TJ 
Localization is described by three terms L m (x), L r p (x) and L r c (x). 

* (0)= E lM°> C i>, 



L m (x) = r_(x) 



E^ 



j ei 



U(x)=T + {x) 



E ~ ^) 



L;(a:,t)=2^Re 



r x (x,t) [ £>i E^-^) 

J'6K K / \ieK„ 



^|c| 2 (cos0±|c|) 2 f x^^H 2 

r ±0) = ^2 ^ *o(a0 + (i - o (x)) ^— 



x ^ -<5 (x) 1 



t+1 



|e| 2 ( cos 2 ^ - |c| 2 ) 



+ (1 - S (x)) 



1 - 



, a x y v/x^r ; 

Then we have the following theorem. 

Theorem 1 (Localization) For k > 1, x G Z + U {0}, r € K K; 
P(X t , r = x) ~ 

i + (— iy+ x 

{i[-l,|c|) (C0S<£) L m (x) +I(_| C |,1] (COS0) L£(x) + I(-| c |,|c|) (cos0) L r c (x,t)} , 

where f(t) ~ means f(t)/g(t) — > 1 (i — > oo). 

We see that in many cases the quantum walk on J K exhibits localization. Localization does not occur only 
in the following two cases, "EjW'j — Vv) = for any r and coscf) > |c|" and "J^. -0j = and cos0 < — |c|". 
Only the symmetric initial state (i.e., ipi = for any «) satisfies the first condition. Moreover L r c (x,t) is 
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Figure 2: Comparison between theoretical values with circles and numerical estimations with crosses of 
W t .3 with C = e^H. The probability P(X tfi = 1) is plotted. The initial state is ip = e ll07r / 180 /V3, 
ip! = e a ^/ im /Vi, ip 2 = e 134077 / 180 /^ and (a) (p = 0,(b)<p = 40tt/180, (c) ip = 50tt/180, (d) tp = 80tt/180. 
Since |c| = l/y/2, ip = 45tt/180 is a critical point for the oscillatory behavior. In the large figure of (c), 
theoretical value arc omitted. 



an oscillatory term, so the probability oscillates if L r c {x,t) exists. The probability P(X tt o = 1) is shown in 
Fig. 2, where we choose the local coin operator as e lip H. From Theorem 1, the condition for the existence 
L^(x,t) is — |c| < cos(p < \c\. Therefore, in this case, the oscillation emerges when 7r/4 < <p < 3tt/A. Remark 
that from Theorem 1 we can see the following relation, 



r x (x,t) E^ E £(V>i-iM 



u'en 



0. 



This means that the oscillation disappears when we take the probability summed over all vertices with a 
same distance from the origin. In addition, since P(W t . K = 0) = Ylj P(Xt,j = 0), the probability of the 
origin does not oscillate for any condition. We also find that the distribution has an exponentially decay 
with x from Theorem 1. The probability -P(^ioooo,o = %) is shown in Fig. 3. 

In order to state the weak convergence theorem, at first we define some parameters depending on the 
initial state. For r € K K , we put 



W) = \lpr 



E 
E 



(3.5) 
(3.6) 

(3.7) 



j,fceK K \{T-} 



Next, we introduce the following notations. Terms C m and C r d are delta measures which are caused by local- 
ization and Cj(x) is a weight on density function /k(x) which is formed a typical shape of one dimensional 
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Figure 3: Comparison between theoretical values with circles and numerical estimations with crosses of 
W t<3 with C = e iv H. The probability P(X 100QQfi = x) is plotted. The initial state is Vo = e llchr / 180 /V3, 
ip! = e i307T / 180 /%/3, ip 2 = e 1340 ^/ 180 /^ and (a) ip = 0, (b) ip = 40tt/180, (c) <p = 50tt/180, (d) p = 80tt/180. 



quantum walks. 



_ b£|c|(|c| - cos< 



C r d {x) 



2K_ 



jeK K 



_ b 2 K \c\(\c\ +cos( 
2K~ + 



ri^W + 2Re(r 2 (x)g 2 r (v>)) + r 3 (x)gg(t/>) 

- (1 - x 2 ) sin 2 - (1 - x 2 ) sin 2 <t>) ' 



T\{x) = 4o K |c|(|a| 2 - x 2 ) cos</>sin 2 <j> 

+ (a 2 + 2a K \c\ cos <p + l) (l + |c| 2 - 2|c| 2 cos 2 - (1 - x 2 ) sin 2 
r 2 (a;) = -2o K |c|(|a| 2 -x 2 )ie^cos0sin0 

+6 K (a K + |c|e'*)(l + |c| 2 - 2|c| 2 cos 2 - (1 - x 2 ) sin 2 0), 
r 3 (a;) = 6 2 (1 + |c| 2 - 2|c| 2 cos 2 <p - (1 - a; 2 ) sin 2 0). 



The weak convergence theorem is derived as follows. 

Theorem 2 (Weak convergence) For k > 1, re K K , as t — > 00 

TTie Zimzt measure is defined by 

Pw( x ) = {l[-i,\c\)(cos4>)Cm + I(-\c\,i] (cos<p)C^}S (x) + C r d (x)f K (x), 

where 



f K {x) = 



I[0,a)(x)y/1 - M : 



ir(l-x 2 )y/\a\ 2 -x 2 
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Figure 4: Comparison between theoretical values with line and numerical estimations with dots of Wt,3 
with C = e lv H. Density function C%(x)fic(x) and scaled numerical values at time 2000 are plotted. In 
theoretical values, the delta measure is omitted. Since localization has an exponentially decay, the scaled 
numerical values corresponding to localization converge to the delta measure at the infinite time. The initial 
state is Vo = e 4l07r / 180 /V3, ipi = e 130 ^ 180 /^, ip 2 = e 1340 ^ 180 /^ and (a) tp = 0, (b) ip = 40tt/180, (c) 
if = 50tt/180, (d) ip = 80tt/180. 



The function C2(x)fic(x) and the scaled numerical values are shown in Fig. 21 Same as other one-dimensional 
cases [13H17] , this distribution has scaling order t and the typical density function of quantum walks Jk (x) . 
The delta measures C m and C r v are caused by localization, i.e., C m — J2 X L m (x) and C r v = J2 X L r p {x). 

The above expressions of both theorems seem to be complicated, however for the following cases, they 
are written in simpler forms. 



Corollary 3 {<p = 0) For x e Z+ U {0} ; r G IK K and <f) = 0, 



P(X t , r = x) 



1+ c 



6 (x) + (1 - So(x)) 



1- Id 



l + \c\J \ I - M 

and X t , r /t converges weakly to a limit measure pyp Q as t — » oo, where 



2(i + |c[ 



5o(x) + 



\c\ - 1)^V + b K J2 $i 

3&.n 



\a\ 2 \A 



x 2 /k(x). 



Proof : From Theorem [TJ when 4> = 0, we have P(X t , r = x) ~ (1 + (-l)* +:E )i^(a;)/2. Also Theorem [2] 
implies r^x) = |a| 2 (a 2 + 2a K \c\ + 1), r 2 (x) = \a\ 2 b K (a K + \c\),T 3 (x) = |a| 2 6 2 , and K + K_ = |a| 4 . After some 
calculations of 0\{ip), (^(VO) 0^{xp), we have the desired conclusion. □ 
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From the definition of Wt, K j Wt i is simply a quantum walk on a half line with a reflecting wall on the 
origin. In the next corollary, we denote X tj o = X t to assert that the walk is defined on a half line. 

Corollary 4 (Half line) For x£Z + U {0}, 

'1 + (-l) t+x 



P{X t = x) ~ I[-l,\c\) (cos 



2 



4|c| 2 (cos0-|c|) 2 | /|ar ' ' 

do (a;) + (1 - d (a;)J — 



if 2 uv ' v ,n " UJ V /£"_ 



and Xt/t converges weakly to a limit measure px as t —¥ oo, where 

t\ T I ^ 2|c|(|c|-cos 0) ^ 2(l-|c|cos^) 2 

M*) = J [ _ 1 ,| C]) (COS0) — S (X) + K __ {1 _ x 2 )sin 2^ fK(x). 

Proof: For k = 1, we have &i = 2/« = 2, oi = h - 1 = 1, 0?(t/>) = 1 and 6>§(V) = ^(VO = 0. From Theorem 
[TJ we get Lp(x) = and L®(x, t) — 0, thus we should consider only L^ n {x) as localization factor of X t . From 
Theorem we have (7° = and 

ri(ar) = 2(1 - |c|cos</>) {X+ - (1 - x 2 ) sin 2 0} . (3.8) 
Combining C m and C^(x) with Eq. (|3.8[) implies px(x). □ 

Remark that we get another proof of Corollary!?] by considering Wt. K with the symmetric initial state. 

We can adopt Corollary [3] for V t K with no perturbation, i.e., 5=1. In <f) — case, the formula for the 
rth half line is directly expressed by V'o, 0i> • • • > V^-i instead of 9\{tp) , O^i^) , O^ijj) . Both cases of = 
and half line, the oscillatory term _L£(iE,i) appearing in Theorem [T] vanishes. 



4 Proofs of Theorems [T] and [2] 

In order to prove Theorems [1] and [5] we consider a reduction of Wt, K on a half line. For W ttK with arbitrary 
initial states, we can not construct the reduction of the walk directly, since the states with the same distance 
from the origin have different amplitudes. To solve this problem, we introduce W' t K which is a quantum 
walk with an enlarged basis of Wt. K - After that, we construct X| as a reduction of W[ K on a half line. To 
analyze X£, we give the generating function of the states. By using it, we obtain the limit states and the 
characteristic function of Wt, K - 



4.1 Reduction to a half line 

Let ^t{x) be the state of the quantum walk Wt, K at time t and position x. We denote the initial state 'Po(O) 
as if> = X^eK„ V'jlO, £j). Now we rewrite if> using a new orthogonal basis {\e'j);j £ K K } as 

where Tvk is the identity operator on Jif^ K > and we defined as 

A(0) = ^r{e' r \®Iw- 

Now let M" be a Hilbert space spanned by an orthonormal basis { | e ^- ) ; j £ K K }. Then we define W[ K 
as a quantum walk on 3%" <S> Jif^ K > with the evolution operator Uj — (I K ® Sj)(I K <8> i*j) — I K ®Uj and 
the initial state 2j£K« l e j)|0> e j)' where Ik is the identity operator on M 7 ' . Let Iq = {eo, ei, • • • , e K -i} and 
l-x — {Up, Down} for a; £ V(I K ) \ {0}. Then the state of quantum walk W[ K at time t and position x is 
written as ^' t (x) = X^gK uei a t( e j' x i u )\ e 'j) \ x > u )> where a' t (a,b,c) is the amplitude of the base \a,b,c) at 
time i. From the construction we obtain the following lemma. 
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Lemma 2 (Enlarging basis) For any t > and x € V(J K ), 

* t (a:)=A(V»)* / t (x). 

Proof : We show the equation by induction with respect to t. At t = 0, by definition of A(t/>), it is trivial. 
For fixed t > 1, we assume # t (x) = A(^>)*f(x), then for a; G V(J K ), 

ps*^) - E «t«-.*.«)iej>i*.«> 

= UjA(i/>) 22 a((e^,a;,w)|ej-)|x,u) 

= X V'j-^( e j> a; > U )l a; J ' U ) 

= X ^c4(4,x,7x)t7j|a;,u) 

= A(V>) E a^x,^}^^}) 

= A^)^(^)- 

This relation holds for any x, so we conclude \l/ t _|_i(x) = A('0)^'j , 1 (x). □ 

For x e V(J«), we define the probability of = x" by P(W/ !fi = x) = \\k(ip)%(x)\\ 2 . Then it follows 

form Lemma[2]that P(W[ K = x) = P{W t , K = %) for any x € V(I K ). For W[ K , the information of the initial 
state is covered by A(i/>). In other words, for any initial state of Wt, K , h is enough to consider the initial 
state X^gk« l e j)|0' e j) 011 W/,k- Consequently, the states of the quantum walk W t K have a good symmetry, 
so we can treat the reduction of the walk. 

Now we introduce X* as a reduction of W' tK on a half line. Here X% is defined on a Hilbert space 
generated by the following new basis. For all I £ {Up, Down} and x € Z + , 

\Own,0,e) = £ I4-.0.6,-), 

| O</ier; 0, e ) = — L= ^ J] |Cfc,0, Ci ), 

|Ou;n,x,0 = X ki'^'W'O) 

\Other,x,l) = J— ^ X l<4A(XM)- 

On this basis, we obtain the one-step time evolution. For x 6 Z + , 

|O«m,0, e) — > a K \Own, 1, Down) + ■y/K— Xb K \ Other, 1 , Down) , 
\Other, 0, e) — > \//t — l& re |Ou>n., 1, Down) — a K \Other,l, Down), 
jj, \Own,x,Up) a\Own,x — l,Up) + c\Own,x + 1, Down) , 
J ' \Other,x,Up) — > a\Other,x — I, Up) + c\Other,x + l,Down), 
\Own, x, Down) — > b\Own, x — 1, £/f>) + d|0wn, x + 1, Down), 
\Other, x, Down) —> b\Other, x — l,Up) + d\Other, x + 1, Down). 
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The subspace generated by this basis is invariant under the operation U'j . Moreover the initial state of W[ K 
can be written as \Own, 0,e). Therefore we can write the evolution operator of A t * as U* H = F^S* H . The 
coin operator is defined by 



F* H = 



^/k— lb K 



\fn-lb K -a K 



<g> |0> <0| ® 1 + /2 ® ® C '- 



For m € {Own, Other}, I <E {Up, Down}, the shift operator S 1 * is defined by 



Sff\m, 0, e) = |m, 1, Down), 

\m, 0, e), Z = f/p, 

|m, 2, Down), I = Down, 



S* H \m,l,l) 



Sh\m,x,l) = \ K^-J'^)' ; = x > 2 . 
I |m, x + 1, Down), t = Down, 

Throughout this paper, we put |Oum) = T [l,0] and \Other) = T [0, 1]. An expression of A t * using weights is 
shown by Fig. [5] and Eqs. (|.18[l - (|.20jl in Appendix. 

Let ^1(x) be the state of the quantum walk X% . Now we define for r € K K , 



A r ('0) = ?/v(Own| + Vk - 1 ^2 ipj(Other\ <8> J 

j£K K \{r} 



Then we introduce A t>r whose probability of u X t , r = x" is defined by 

P(X ttr = x) = \\A r (iP)y* t (x)\\ 2 . (4.9) 
This probability is described by the state of Wt K in the following, 



P(X t>r = ar) 



|a t (0,e r )| 2 , i = 0, 
||^I't(/i r .(a;))|| 2 , otherwise. 



Hence the relation between the probabilities of "Wt, K = h r (x)" and "X t . r = x" is obtained as 

f E P(Xt tj = o), x = o, 

P(W t , K =h r (x)) = i ieK, 

[ P{X t>r — x), otherwise. 

Note that £ j6Kk Exe{o}uz+ =«) = !• 

In Subsections 14.21 and 14.31 we analyze ^ (x) by the generating function. 

4.2 Proof of Theorem Q] 

We compute the limit state of A t * from the generating function which is defined by 

oo 

**(x; z) = V" ** (a;)^' = V" a*(Z,a;,m;2;)|Z,x,m). 

t=0 le{Own,Other} 
rn£{Up,Down} 
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From Appendix, we see that there exists < r± < 1 so that for any z with \z\ <r%, 

a* (Own x Utrz) - I ~ltMz) ~ az){u{z) + a K )<5>{x; z), x > 0, . , 

a [Uwn, x, Up, z) | _ {fi{z) + anMzMx . g)j x = 0i ^.iu; 

a*{Other,x,Up;z) = { -^iK^)-az)b^{x ]Z ), x > 0, (4 n) 
^ — vk — lb K n(z)9{x; z), x = 0, 

~*/^ \ ( —z(ri(z) + a K )^(x;z), x > 0, ^-,^\ 

a*(Owti, x, Down;z) = I q a; — ^ ^ 



n \ / ~z\/k- lb K $(x;z), x>0, fAti\ 

a (Other, x, Down; z) = < _n (4-13) 



( a J 4(1- C 2)(z2_ w ^)( z 2_ w 2 ) > ( 4 - 14 ) 



where 



w , Az 2 + 1 - v /A 2 z 4 + 2A(l-2|a| 2 )z 2 + l 

A(Z) = 2dz ■ 

, , d\(z) - Az 

H(z) = z, 

c 

v{z) = (l + Az 2 ) 2 -4A|a| 2 z 2 , 
r)±(z) = 2c ± 1 T Az 2 , 

2 _ C(l ± C) 

W± - T A(|a| 2 -l T c)' 
A = ad — 6c. 



Note that | = (1 ± c)/(l ± c) =1. From Cauchy's theorem, we have for < r < ri < 1, 

Therefore as t — > oo 

- V*(x) ~ Res(^*{x;z),w + )w+ {t+1) + Res(**(a;; z), -w+){-w+y {t+1) 

+ Res(\&*(x; z), w_)wZ (t+1) + Res(**(a;; z), , 

where Res(/(z), w) is the residue of /(z) for z = w. Taking the residues of the generating function, we can 
compute &*(x). After some calculations with Eq.(9) and &*(x), the proof of Theorem 1 is complete. 

4.3 Proof of Theorem H 

In order to prove Theorem we calculate the Fourier transform of the generating function as ^*(s;z) = 
J2 X ^* ( x 'i z)e lsx by Eqs. (|4.10p - (14.13[) . Then we obtain the characteristic function from the following relation 



f27T 



2tt 







(ArW)yUx),KWW(y)y k{x - v) e i£x 



dk 

27 



\A r W)4>* t (s),ArWW(s + 0)^, (4-15) 
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where (u, v) is the inner product of vectors u and v. 

Now we write the Fourier transform of the generating function as 



ty*(s;z)= a* (l,m; s; z)\l,m) . 

i e 

m i 

From Eqs. (pnD|) - (|i33"]) . we have \?*(s; z) as 



l e {Own, Other} 
m £ {Up, Down} 



where 



a* (Own, Up; s; z) = \ —fJ,(z) + — (\(z) - az)$ 2 (s; z) J (a« + ^(z))<f>i(s; z), 
a* (Own, Down; s; z) = — H (A(^) — az)$2(s; z) ) %/k — l& re $i(s; z), 

V ac y 

a* (Other, Up; s; z) = z(a K + jj,(z))3>\(s; z)$2(s; z), 
a* (Other, Down; s; z) = zyn — lfe K $i(s; z)$2(s; 2), 

wiw 2 0) + v / K^)) (f- ( z ) - V^)) 

$1 (s- zl = ~ - — - - 

U ' ^ 4(1 -C 2 ) (Z 2 - W 2.)(Z 2 - W 2) 



(C( S ;^)-V^)) 



e lfe 



2A(z - v + (s))[z - v-(s))' 

C(s;z) = 2ae~ ls z~l - Az 2 , 

. , ae" ls + aAe ls ± J(ae- ts + aAe ls ) 2 - 4A 
*>±(s) = ' 



2A 

Here we can rewrite v±(s) as 



(s) = e~ lp (\a\ cos 7 (s) ± ^M 2 cos 2 7 (s) - 1) = e - ip e ±i6 ^, 



where we take A = e 2lp , a — \a\e la , j(s) = s — a + p, cos9(s) = |a|cos7(s). Note that |u±(s)| = 1. Now 
||**(s;z)|| 2 < oo for < \z\ < n, we can rewrite *S>*(s;z) = J2t>o *t ( s ) 2 '- So we nave for < r < ri 

If ~ , , dz 



Therefore we get the Fourier transform of the state ^ (s) as follows: 

- ~ ^ + ( s )( w+ )- (t+1) + i>- w+ ( s )(-w + y^ 

+ i, w _ ( s )( w _r( t+i > + ( s )(-w^y^ 

+ i; v+ ( S )(v + (s))-^ +ip v _(s)(v-(s))-( t+1 \ (4.16) 

where ip± w ±(s) = Res(\&(s; z); ±w±) and Vv ± (s) = Res(^(s; z); v±(s)). 

Finally we compute the characteristic function by Eas. (|4.15p and (|4.16D . Now we have 

27r dn 

(\\A r (^ w+ ( s )\\ 2 + \\A r (^ w+ (s)\\ 2 ) — = c;, 

27r dn 

(\\A r m w _(s)f + \\A r (^ w _( s )\\ 2 ) — = C m , 

2 * ds 
((A r (^)ip w+ (s),A r (^- w+ (s)) + (A r (^- w+ (s),A r (ij})ip w+ (s))) — = 0, 

27r ds 
((A r (tp)ip w _ (s), A r (tp)i(j- W _ («)) + (A r (ip)tp- W _ (s), A r (tl))tp w _ (*))) — = 0. 
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Noting that e i(t+VH*+t/t) = e i(t+i)e(s)+i£h(s)+o(t 1 ) ) w h er e /i(s) = dO{s)/ds, the above equation and Ea. (l4~T5l) 
with the Riemann-Lebesgue lemma imply 



lim E 

t— ¥00 



oitXt, r /t 



2tt 



= c; + c m + I e-^pts)— 







2tt 



(4.17) 



where p(s) = ||A r (^>)^ 0+ (s)|| 2 and q(s) = ||A r (t/>)^i„_ (s)|| 2 . Moreover, from a change of variable for last two 
terms in Eq. (|4.17[) . we have 

f e^ h(s) p{s) + e l ^ s) q(s)) — = [ e* x w(x)f K {x)dx. 
After some calculations for w(x) with p(s) and q(s), we have the desired conclusion. 



5 Summary and discussions 

We introduced a quantum walk with an enlarged basis to consider a reduction of quantum walks with 
arbitrary initial state. This method is based on an idea canceling the asymmetry caused from initial state 
by a new tensor product. From our results in this paper, we discuss two interesting points. First, we found 
the oscillating probability as localization. From Theorem [TJ the oscillatory term is expressed by L r c (x,t). 
We can see that this term vanishes with some initial states or local coins. For example, we consider Vt jK , 
which is a quantum walk on T K with local coin operators G K with additional complex phase c at the origin. 
If c = 1 Corollary [3] implies that the oscillation does not occur with arbitrary initial state. Also if the initial 
state is symmetric, the walk is reduced on a half line. Then it follows from Corollary H that no oscillatory 
behavior arise with any complex phase c. Thus the initial state and differences on complex phase of local 
coins are important factors for the oscillatory behavior on localization. Especially in quantum walks on 
the one-dimensional lattice with homogeneous local coins, localization does not occur [TSl - fTB] . If localization 
occurs with perturbations of local coin operators on the one-dimensional lattice, there seems to be a condition 
that an oscillatory behavior arises in localization. Second, Wt. K has the scaling order t and the limit measure 
has the density function Jk{x) which is a half- line version of one appearing in the quantum walk on a 
line [T3HT5] . This is a typical property of quantum walks [8l[T6j[T7]. To show the universality of the limit 
theorems for quantum walks is one of the interesting future's problems. 

Acknowledgments. We thank Noriko Saitoh and Jun Kodama for useful discussions. N.K. is supported 
by the Grant-in-Aid for Scientific Research (C) (No. 21540118). 
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Appendix 








Q P Q 

Figure 5: Quantum walk with enlarged bases X£ on a half line 

We calculate the generating function of ^*(x) by using the method in [53]. To simplify notations, for 
I G {Own, Other}, we denote | Z , , e) = | Z , , Up) and construct Down) as a dummy base, which always 
has value as its amplitude, so that the local coin operator on the origin has 2x2 matrix. To indicate the 
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X 



X\ 

1 













/ V \/Q 



t 



;(0^0;n) ri+T2 + l 

Figure 6: 3(0 ->■ X\\t\ + t 2 + 1) 



evolution operator of the walk, we use an expression using weights (see Figj5j) , where 

a K 



Vk- lb K -a K 



a b 






1 

" 
c d 



' 1 " 




' 1 " 











(.18) 
(.19) 
(.20) 



Q 

p = / 2 ® 

We define the generating function for the state by 

oo 



In order to compute ^*(x; z), we first define the transition amplitude 3(0 — > x; r) as the weight of all paths 
starting from ending at x after r steps, and 5(0 — > x; r) as the weight of all paths on another walk defined by 
Q' = Q. For example, 5(0 -> 2; 4) = QPQQ+PQQQ+QQPQ and 3(0 -> 2; 4) = QPQQ+PQQQ+QQPQ. 
From 3(0 — > 0; n) and 3(0 — > x — 1; t 2 ), we can obtain 3(0 — > x; t\ + t 2 + 1) as Fig[5] Then we get ^(x; z) 
from the generating function for 3(0 — > x; r). 

We now calculate the generating function for 3(0 — > x;t). Since the first operator should be Q on the 
half line, the weights of paths form Q ■ ■ ■ Q or P- ■ ■ Q. So we express 3(0 — > x; r) as a linear combination of 
Q and R: 



where R = 7 2 



c d 


5(0 — > x; t) is defined by 



3(0 -> x- t) = b q (0 -> x- t)Q + b r {0 -> s; t)R + S (x)S (t)I 2 ® J 2 

and we define & 9 (0 — > x; 0) = 6 r (0 -+ x; 0) = 0. The generating function for 



^ 5(0 -+ x; r)z T = B 9 (0 ->■ x; z)Q + B r (0 -t x; + 6 (x)I 2 ® 7 2 , 

with £9(0 -> x;z) = E^Lo 69 (° -> z;r)z r and £ r (0 -> x; z) = £^ =0 & r (0 ~> x;t)z t . Since the left-hand 
tensor product of P and Q is 7 2 , the generating function for 5(0 — > x; r) corresponds to the result in [23], 
i.e., for sufficiently small z, 

B q (0-^x;z) = (-\(z)\ i x>l, 
a a 



S«(0 -> 0;z) = 0, 
B r (0 -+ x-z) 

X(z) 



d . , A X(z) — az 

-X(z) } — , x > 0, 

a acz 



Az 2 + 1- v /A 2 z 4 + 2A(l-2|a| 2 )z 2 + l 



2dz 



(.21) 
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Here we take A(z) for the smaller solution of the absolute value of 



A 2 W 



Az 



A(z) + -=0. 



(.22) 



Note that for sufficiently small z we can write X(z) by Eq. (j.21l) . Moreover since \a/d\ = 1, we can take 
ro < 1 such that |A(z)| < 1 for \z\ < Tq, Next we calculate the generating function for 5(0 — > 0; r). To do so, 
we introduce a new notation 5(0 — > 0; t; n) as the weight of all paths starting from the origin reaching the 
origin n times before ending at the origin at time r. Now we consider S(0 — > 0; t; 0). For r > 2, we obtain 
5(0^ 0;t;0) as (see Fig0) 



where R 



a K Vk - lb K 
VK-lb K -a K 





(1- 


6 2 {r))P{b r {0^0: t -2)R}Q 


F <J 2 (r)PQ 




{(1- 


-s 2 


(r))a66 r (0-^0;r-2) + 5 2 


(r)6}4 


® 


" 1 



" 



and for r < 2 we define E 


(0^0;r;0) 


^0; 


r;0) 


as 







^ 5(0 -> 0; t; 0)z r = (adB r (0 0; z) + &)z 2 ,R = B f (0 -> 0; z; 0)A 



T = 



Similarly, for r > 4 we have S(0 — > 0; r; 1) as 

5(0^0;r;l)= ^ {(1 - 5 2 ( n ))a66 r (0 0;n) + fe(n)&}£ 

Tl+T2+4=T 

x {(1 - 6 2 (T 2 ))aWF(0 -> 0; r 2 ) + <5 2 (r 2 )&}i?, 
and for t < 4 we define 5(0 — > 0; r; 1) = 0. Thus the generating function for 5(0 — > 0; r; 1) is obtained by 

oo 

H(0 -> 0; r; l)z T = {(adB r (0 -> 0; z) + fo)z 2 } 2 i?/ = (0 -> 0; z; 



T=0 



where i?/ = I 2 



1 




Recursively we have the following formulae: for n > 0, 
l + (-l) 



£ r (0 -> 0;z;n) = 
B fl (0^0;2;ri) = 



{(adB r (0^0;z) + 6)z i } 



2in+l 



1 + (_1)»+1 



{(adB r (0^0;z) + 6)z 2 } 



2l Tl+1 



(.23) 
(.24) 



From Eqs. (|.23p and (l.24[) . we get the generation function for 5(0 — > 0; r) by summing over n. Here 
(dB r (0 — > 0; z) — Az)z 2 = (<iA(z) — Az)z/c, so we see that for z with |z| < ri = min(|c|,ro), 



\{dX(z) - Az)z/c\ 2 < (\dX(z)\ 2 + \z\ 2 )\z/c\ 2 < \d\ 2 + \c\ 2 = 1. 
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Therefore for z such that \z\ < ri, 

oo 

2(0 -> 0; r)z T = B f (0 -> 0; z)i? + B f '(0 -> 0; + I 2 ® J 2 , 

of/r, n \ V^nf/n n \ (dA(z) - Az)z/c 

B (0 -> 0; z) = g (0 -> 0; z; n) = ± {g(rfA(z) - ^ } , 



B fl (0 -> 0; z) = B? ' (° ^ 0; z; n) = y 



n=0 



{(d\(z) - Az)z/c} 
-{(d\(z)- Az)z/c} 



2 

~2 ' 



For £ > 1, 5(0 — > x; r) is written by c(0 — »• 0; r) and 5(0 — » x; r) (see FigE]) as 

S(0^x;r)= S(0->x-l;7a)QS(0->0;Ti)+ib(T)*o(a;)Ja® J a . 



T 1 +T 2 + l = T 



From the generating function for 5(0 — ► 0; r) and 5(0 — > a;; r), we can compute the generating function for 
5(0 — > x; r) as follows: for x > 1, 



]T 5(0 x; r)z T = {B 9 (0 x - 1; z)Q + B r {0 -+x-l;z)R+ 8x(x)I 2 ® / 2 }Q2 

T=0 

x{B f '(0 -> 0; z)i? + £ f '(0 -> 0; z)i? 7 + I 2 ® I 2 } + 5 {x)I 2 ® 7 2 . 
Then we obtain the generating function vp* (x; z) as follows: 

oo 

§*(x;2) = ^S(0^x;t)z^(x). 
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